Quantum size effects on the perpendicular upper critical field in ultra-thin lead films 
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We report the thickness-dependent (in terms of atomic layers) oscillation behavior of the perpen- 
dicular upper critical field H C 2± in the ultra-thin lead films at the reduced temperature (t = T/T c ). 
Distinct oscillations of the normal-state resistivity as a function of film thickness have also been ob- 
served. Compared with the T c oscillation, the H C 2± shows a considerable large oscillation amplitude 
and a ir phase shift. The oscillatory mean free path caused by quantum size effect plays a role in 
Hc2± oscillation. 
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There is a long history of scientific research on super- 
conducting thin films. In particular, theoretical and ex- 
perimental studies have been carried out to understand 
how the film thickness affects the superconducting prop- 
erties. It seems that the reported experimental results 
of thin films can be explained by the existing theories 
of superconductivity 111 |^ yl 14 S la. Q • However, most 
previously studied superconducting films were still rel- 
atively thick, normally over several tens of nanometers, 
and the film morphology was usually poor. If the film 
surface is atomically uniform and the thickness is further 
reduced to several nanometers so that the quantum size 
effects become apparent, a natural question arises: will 
some unexpected new phenomena emerge? In particular, 
does the conventional theory still work? 

Previous theoretical works have predicted many pos- 
sible prominent physical properties modulated by quan- 
tum size effects: electronic structure, critical tempera- 
ture, electron-phonon interaction, resistivity, Hall con- 
ductivity and so on H, II H El E 111 There 
are also some related important experimental results 
BEHmEHIIllIll, such as the T c oscillations 
in ultra-thin Pb films, which are caused by the density 
of states oscillations in confined quantum well structures 
[iil E^j l and by the electron-electron interaction medi- 
ated by quantized confined phonons [ill However, 
the properties of the upper critical field affected by the 
quantum size effect have not been reported in previous 
work. In this Letter, we report our experimental observa- 
tion of the oscillatory H c2 ± through magneto-transport 
measurement of ultra-thin Pb films. The oscillations are 
similar to those of T c but the motivations are more com- 
plex. Besides the factors for T c oscillation, we interpret 
this unexpected phenomena by the oscillatory mean free 
path in ultra-thin superconducting films caused by the 
quantum size effect. 

The 3mmxl0mm sized Si(lll) wafers were used as 
substrates and prepared by the standard cleaning pro- 
cedure to obtain the clean Si(lll)-7 x 7 surface. The 
base pressure of the UHV-MBE-STM-ARPES (Angle 



Resolved Photoemission Spectroscopy) combined system 
we used was about 5 x 10~ n Torr. The Si substrate 
was cooled down to 145K during the MBE layer-by-layer 
growth of the Pb films. The growth rate was controlled 
at 0.2 ML/Min (Monolayer/Minute) and a RHEED (Re- 
flection High Energy Electron Diffraction) was used for 
real time monitoring of the growth. After deposition, the 
sample was warmed up slowly to room temperature and 
transferred to the analysis chamber where the STM and 
ARPES were used to investigate the surface topography 
and the electronic structures, respectively [24)]. For ex- 
situ magneto-transport measurements, all the Pb films 
were covered with a Au protection layer of 4ML before 
being taken out of the UHV system. 

The R — H measurements were carried out in a very 
short time after the samples were taken out of the vac- 
uum. The applied field was perpendicular to the sample 
surface and the temperatures were set near and below T c . 
To avoid trapping flux in, the magnet was discharged to 
zero in oscillate mode and the sample was warmed up 
to 8K before the R — H measurement for each temper- 
ature. Then the perpendicular upper critical field -ff C 2_i_ 
at different temperatures was obtained from the R — H 
measurements at the field where the resistance reached 
half of the normal state resistance Rn- The resistance 
approaches Rn very gradually because of the magnetore- 
sistance effect. So we took Rn as the resistance where 
the resistance variation ratio is within 0.1%. 

Figure 1 shows the R— H curves of a 21 ML sample at 
different temperatures. The arrow points out the defined 
perpendicular upper critical field H C 2±_ at 4.7K. The inset 
of Fig. 1 shows H C 2± vs. temperature for the 21 ML film. 
It shows a perfect linear dependence on T near T c , which 
is a typical property of a superconductor with a high 
value of the Ginzburg-Landau parameter k. The inset 
of Fig.l can be used to determine the zero field critical 
temperature T c by extrapolating the plot to H C 2±_ = 0. T c 
determined in this way is shown as a function of thickness 
in Fig. 4(a). Normally, a direct way of determining critical 
temperature is through the R — T measurement at zero 
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FIG. 1: R-H curve of the 21 ML sample. The magnetic field 
is perpendicular to the sample surface. The black arrow in- 
dicates the determined upper critical field at 4.70K. The in- 
set shows the H C 2± as a function of T for this sample. The 
plot is linearly extrapolated with dashed lines to both high 
and low temperature sides. The measurements were carried 
out with a Quantum Design Magnetic Property Measurement 
System(MPMS-5). 
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FIG. 2: The reduced resistances of Pb films as a function of 
temperature are shown in (a). The resistances are normalized 
by the normal state resistance at T = 8K. (b) shows an 
oscillation of normal state resistivity at 8K as a function of 
film thickness. 



field. We find that the critical temperatures determined 
by both methods show a consistent oscillation behavior 
and the values are quite close for every thickness. 

The reduced R — T curves of Pb films from 21 ML to 
28 ML are shown in Fig. 2(a). The normal state resistiv- 
ity p n oscillation with film thickness at T = 8K is shown 
in Fig. 2(b). The similar resistivity oscillations caused by 
quantum size effect have been reported in single crys- 
talline Pb and Pb-In thin films at T = 1WK 16]. But 
in polycrystalline films, oscillations of the normal state 
resistivity have not been observed although T c has been 
found to oscillate with film thickness [la]. In our ex- 
periment, both T c and p n oscillations are observed. It 
indicates that the quantum size effects show up in both 
superconducting state and normal state but the intensi- 



ties and mechanisms may vary in different way depending 
on sample conditions. 

Figure 3(a) shows H C 2±. as a function of the reduced 
temperature t = T/T c . For every thickness i? C 2J- shows 
a good linear dependence on t near t = 1. H C 2± ver- 
sus film thickness for t=0.90 and 0.95 are shown respec- 
tively in Fig. 3(b). It is shown that with the film thickness 
variation, H c2 ± exhibits an oscillation behavior, which is 
similar to the reported T c oscillation |22|. However, the 
oscillation of H C 2± are ir out of phase to that of T C) i.e., 
peaks appear in the odd layer samples where dips appear 
in the even layer samples, which is opposite to the T c 
oscillation shown in Fig. 4(a). 

In the early theories proposed to understand the mag- 
netic properties of thin film superconductors, the TGS 
(Tinkham, de Gennes and Saint James) theory 00 was 
validated as showing a good agreement with the former 
experimental results 00. According to TGS theory, the 
upper critical fields H C 2± near T c should monotonically 
increase when the film thickness decreases, which can be 
described in the following form 0: 



H c± (T,d) = V2K(T,oo)H c (T)(l + b/d), 



(1) 



where k(T,oo) = 2v / 27r# c (T)A^ o (T)/$ and b = 
3A|(T)Co/8A^(T). Here H C (T) is the thermodynamic 
critical field, Al is the London penetration depth, Aqo 
is the bulk weak field penetration depth, <I>o is the flux 
quantum ($ = hc/2e = 2.07 x l()- 15 Wb) and d is the 
film thickness. In Fig. 3(b), the dashed lines, calculated 
using Eq.(l) and the related parameters in previous work 
with film thicknesses appropriate to our samples, show 
the same tendency as the experimental curves if the os- 
cillations are ignored. The measured H C 2± values of our 
samples are about three times larger than the calcu- 
lated values (note the different scales on the two sides 
of Fig. 3(b)), which may be caused by stronger interface 
or impurity scattering in our films that gives rise to a 
large resistivity, thus large H C 2±_ (see discussion below). 
The linear dependence on t shown in Fig. 3 (a) also gives 
an information that for a given film thickness, the tem- 
perature dependence follows reasonably well with Eq.(l) 
whether that particular film is at the peak or valley of 
the H C 2± oscillation. 

The TGS theory above includes surface scattering ef- 
fects but does not consider the quantum size effects that 
occurs in ultra thin films. The absent H C 2± oscillation 
from TGS theory means that the thickness depended 
quantum size effect is the original source of the H C 2±_ os- 
cillation. According to the G-L (Ginzburg-Landau) the- 
ory, H C 2± is determined by the in-plane coherence length 
In a three dimensional anisotropic superconductor, 
the perpendicular upper critical field near T c is given by 
El 0: H C 2± = TTfij- Our ultra-thin films are thinner 

than lOnm which is much smaller than the Pipard coher- 
ence length of a bulk Pb superconductor (£,Q Ulk — 83nm), 
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FIG. 3: (a) shows the perpendicular upper critical field versus 
the reduced temperature t. The oscillation behavior at t=0.90 
and 0.95 are plotted in (b). The dashed lines correspond to 
the calculated results using Eq.(l). 
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FIG. 4: (a) shows the oscillation behavior of T c with film 
thickness, which is defined by the way shown in the insert of 
Fig.l. The rescaled p n T c variation is shown in (b), which is 
defined in the following way: Ap n T c = (p„T c — p' n T' c ) / p' n T' c , 
where p' n T' c is the value of p n T c for the 21ML film. 



we can use the quasi two dimensional formula [2£ 



dH c2 



d(T/T c ) 



$0 



(2) 



For the linear dependence on t near t = 1 shown in 
Fig. 3 (a), H c2 ± has the same oscillation behavior with 
thickness as that of — ( a*' l ) a t a certain t. The sys- 



tem should be considered as a dirty-limit superconductor 
because of the strong scattering. For dirty superconduc- 
tors near T c , ss £ l where £o is the Pipard coherence 



length and / is the mean free path for a film [3, ■ Ac- 
cording to BCS theory, £ oc 1/T C , therefore we can get 
H c2 ± oc T c /l at a certain t. In Fig. 3(b) and Fig. 4(a), 
it is shown that the oscillation amplitude of H c2 ± and 
T c are about 40% and 10% respectively. On the other 
hand, the mean free path I and the normal state resistiv- 
ity p n have the following relation: I oc l/p n , from which 



we can derive H c2 ± oc p n T c - In Fig. 2(b), the p n oscil- 
lation shows a big amplitude about 60% and the same 
phase as H c2 ±. The rescaled variation of p n T c is shown 
as Ap n T c in Fig. 4(b), which fits well with the oscillation 
behavior of H c2 ±. It implies p n oscillation dominates 
over the T c oscillation in H c2 ± and gives rise to a tt phase 
shift between T c and H c2 ± oscillations. In earlier works, 
the oscillatory conductivity in quantized thin films has 
been observed [16J and the effects of impurity and sur- 
face and interlayer roughness on quantum size effects in 
thin films have been discussed [lZ . [2(| . Though at the 
moment we do not have a complete answer to the oscil- 
lations of p n with thickness for our films with atomically 
uniform surfaces, the previous experiments on the layer- 
spacing oscillation ji^l provides a strong indication that 
the modulation of the interface roughness with thickness 
may play an important role. In that experiment they 
found that the interlayer spacings oscillate with a period 
of quasi double layer and even-monolayer samples have 
shorter interlayer spacings. This is also supported by the 
binding energy modulation observed H2,|^|. It indicates 
the lattice feels at home with the conduction electrons 
for the even-monolayer samples while it is not so for the 
odd-monolayer samples. The unaccommodating lattice 
and conduction elections in the odd-layer samples could 
induce some lattice distortion and therefore enhance the 
interface roughness. This enhanced interface roughness 
must induce a higher resistivity. 

We believe the experimental findings in our ultra-thin 
films are due to a variety of combined quantum size ef- 
fects from ultra-thin film thickness. The quantum size 
effect can show up either as a modulation of the in- 
terface roughness induced by the interlayer spacings as 
well as the modulation of the phonon modes and the 
electron-phonon couplings which both affect the normal 
state transport properties of the samples, of course, also 
causing the wave vector quantization along the thickness 
direction. Under the circumstance, only the components 
of electronic wave vector in the surface plane, i.e., the 
x-y plane, have a continuous distribution. Therefore, 
the electron density distribution is rather inhomogeneous 
along the z direction. The modulation of the electron- 
densities may further feedback to the electron-interface 
and electron-phonon scattering processes and therefore 
to the mean free path. Another relevant issue is that 
the G-L theory is only a mean-field theory in which all 
the short-distance fluctuations are integrated out. For 
our ultra-thin films, to give an adequate description of 
all the electronic states and the scattering processes, we 
must go back to the microscopic theory of BCS super- 
conductivity within the subband framework and derive 
the multi-band G-L theory. The G-L order parameter 
^ perpendicular to the film is limited to quantized val- 
ues and may also show modulation with the interlayer- 
spacing modulation. Each subband may have a different 
value of coherence length £ in the x-y plane, namely t; nj d, 
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where n is the subband index and d is the number of the 
monolayers. In general, H C 2± is determined by a matrix 
equation with m being the size of the matrix in which m 
is the number of subbands below the Fermi energy. In 
the limit that one of the £ nj d is much smaller than all the 
others, H C 2±_ is predominantly determined by this min- 
imum value, which could be much higher than that of 
the bulk. The story here is similar to that of the newly 
discovered superconductor MgB2 , where only two bands 
are involved [27| . If the film becomes thicker, the number 
of subbands will increase. The interaction of subbands 
will weaken the quantum size effects and the coherence 
length will be close to the average one. The oscillation 
behavior of H C 2± will eventually disappear beyond a large 
thickness. 

In conclusion, a large oscillation of H C 2± in the ultra 
thin lead films are observed as a function of film thick- 
ness. The H C 2± oscillation is opposite to that of T c in 
phase and cannot be simply attributed to the modula- 
tion of the density of states and T c . A large value of 
H c2 ± is also observed. Considering the interface and sur- 
face scattering and the modulation of coherence length 
and mean free path induced by the quantum size ef- 
fect, a possible mechanism is proposed to explain both 
the anomalous oscillation of H C 2± and its large value. 
We believe that a quantitative description for the find- 
ings in our experiments must be based on the combined 
quantum size induced modulation effects on the inter- 
layer structures, electronic structures, phonons, electron- 
phonon and electron- interface scattering processes. Fur- 
ther consideration about the flux dynamics is also neces- 
sary by including the interface and surface scattering ef- 
fects and two-dimensional fluctuations in the multi-band 
G-L theory. 
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